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ABSTRACT 


During  the  contract  period  the  main  emphasis  was  on  the  observation 
of  crustal  tilts  associated  with  earthquakes,  their  directions,  amplitudes 
and  propagation  velocities.  Data  from  the  seismology  laboratory’s  three 
borehole  installations  at  GLM,  PAX,  MCK,  the  COL  World  Wide  Standard 
Station  and  the  ALPA  (Alaska  Long  Period  Array)  have  been  analyzed  for 
tilt  steps  associated  with  local  earthquakes.  Tilt  steps  from  distant 
earthquakes  observed  at  GLM,  PAX  and  MCK  and  data  presented  in  the  ixtera- 
ture  are  discussed.  Observed  tilt  steps  cover  a  distance  range  from  8 
to  12600  km  and  a  magnitude  range  from  2  to  8.  Tilt  step  propagation 
velocities  have  been  observed  from  1.2  to  4.0  km/sec,  and  those  corre¬ 
sponding  to  arrival  times  of  teleseismic  S  waves.  Tilt  directions, 
amplitude-  and  velocities  observed  at  several  stations  simultaneously 
for  the  same  earthquake  are  internally  consistent  and  are  likely  to  de¬ 
pend  on  the  tectonic  environment  of  the  observing  station  in  addition 
to  the  focal  mechanism.  The  low  velocity  of  1.2  km/sec  is  consistent 
with  a  plastic  wave  propagation.  Other  velocities  are  close  to  the 
Rayleigh  wave  short-period  group  velocity.  The  tilt  amplitude  for  fixed 
epicenter  distance,  focal  depth  and  stress  system  depends  logarithmically 
on  the  magnitude,  with  a  possible  cut-off,  but  are  too  large  when  com¬ 
pared  to  existing  theoretical  dislocation  models  for  a  single  layered 
homogeneous  isotropic  elastic  crust.  Alternative  solutions  in  terms  of 
more  realistic  crustal  layering  and  including  plastic  deformation  under 
preexisting  tectonic  stress  probably  could  account  for  a  number  of  ob¬ 
servations  . 

The  operation  of  a  number  of  short-period  telemetry  stations  has 
been  discontinued,  new  ones  installed  and  one  was  relocated. 
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A  seismicity  map  was  published  for  the  period  1968-1971,  based  on 
recordings  from  the  telemetry  net  and  the  results  related  to  the  geology 
of  Alaska.  The  map  and  its  geology  related  discussion  proved  very  use¬ 
ful  to  a  number  of  state  agencies  (such  as  the  Highways  Department)  and 
aroused  much  interest  from  private  sources  as  well. 
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TELEMETRY  NETWORK 

A  number  of  changes  in  the  telemetry  network  have  been  made  during 
the  contract  period.  The  operation  of  BIG,  PJD,  SVW  and  BLR  have  been 
discontinued.  MCB  was  replaced  by  ENG,  and  CCB  and  FYU  have  been  newly 
established.  The  following  table  indicates  all  seismic  stations  operated 
by  the  Geophysical  Institute,  though  not  all  under  this  contract  or  by 
the  seismology  laboratory. 


Station  Name 

Augustine  Island 
Broad  Pass 
Birch  Hill 
Clear  Creek  Butte 
Elvey 

Engineer  Hill 
Fort  Yukon 
Gilmore 
Hepp 

Hurricane 

McKinley 

Mercer 

Minitrack 

Nenana 

Paxson 

Remote 

Sheep  Mtn 

Sunshine 

Scotty  Lake 

Tan  ana 


Code 

Latitude  (N) 

Longitude  (W) 

Elevation 

Components 

AGI 

59  22.8 

153  25.2 

580 

SPZ 

BDP 

63  13.87 

149  15.56 

709 

SPZ 

BRH 

64  51.89 

147  38.42 

330 

SPZ 

CCB 

64  38.88 

147  48.32 

183 

SPZ 

ELV 

64  51.62 

147  50.91 

180 

SPZ  Low  Gain 

ENG 

64  44.11 

147  02.76 

220 

SPZ 

FYU 

66  33.96 

145  13.90 

137 

SPZ 

GLM 

64  59.24 

147  23.34 

820 

SPZ,  LP,  Tilt 

HPP 

64  47.43 

147  57.55 

170 

SPZ 

HUR 

62  58.74 

149  38.77 

510 

SPZ 

MCK 

63  43.94 

148  56.1 

610 

SPZ,  LP,  Tilt 

MCR 

63  53.71 

149  Oi.56 

456 

SPZ 

MIK 

64  52.34 

147  49.68 

160 

SPZ 

NEA 

64  34.63 

149  04.63 

365 

SPZ 

PAX 

62  58.25 

145  28.12 

1130 

SPZ,  LP,  Tilt 

RON 

62  41.47 

150  12.22 

470 

SPZ 

SCM 

61  50.00 

147  19.66 

1020 

SPZ 

SSH 

62  10.00 

150  04.66 

100 

SPZ 

SCT 

62  19.15 

150  17.83 

140 

SPZ 

TNN 

65  15.4 

151  54.7 

504 

SPZ 

3 


EARTHQUAKE  EPICENTERS  IN  INTERIOR  ALASKA,  1968-1971 
Contribution  by  L.  Cedney,  L.  Shapiro,  D.  VanWormer  and  F.  Weber* 

On  the  accompanying  nap,  somewhat  over  1500  epicenters  are  repre¬ 
sented.  These  data  were  collected  by  the  Geophysical  Institute’s  U'le- 
meteKed  seismic  net  during  the  four-year  period  1968-1971,  and  the  cal¬ 
culations  giving  the  earthquake  parameters  were  performed  by  the  University 
of  Alaska's  IBM  360-40  computer.  No  distinction  is  made  as  to  earthquake 
magnitude,  although  it  is  likely  that  none  smaller  than  magnitude  2  are 
depicted,  and  that  the  heaviest  concentrations  represent  aftershock  zones 
of  larger  earthquakes.  An  exception  to  the  latter  statement  is  the  dense 
concentration  in  the  Mt.  McKinley  area  (63°N,  151 °W),  where  very  many 
small,  intermediate  depth  earthquakes  occur,  but  earthquakes  of  over 
magnitude  5-1/2  are  rare.  Other  clusters  of  note  occur  at:  64.6°N,  147. 0°W, 
the  area  of  the  Salcha  earthquake  of  July  1937  (m  -  7-1/4);  64.8°N,  147. 4°W, 
the  aftershock  zone  of  the  Fairbanks  earthquakes  of  June  1967  (m  -  6-1/4) ; 
and  65.5#N,  150. 0°W,  the  aftershock  zone  of  the  Rampart  earthquake  of 
October  1968  (m  »  6-1/2) . 

The  correlation  of  epicenters  with  mapped  faults  is  clear  in  the 
case  of  the  Rampart  Fault,  the  Denali  Fault  (extending  east  to  west  across 
the  bottom  of  the  map),  and  several  faults  of  lesser  magnitude.  It  appears 
likely  that  unmapped  faults  may  exist  in  other  areas  where  obvious  epi- 
central  alignment  occurs,  particularly  in  the  Fairbanks  area. 


*  U.  S.  Geological  Survey 
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CRUSTAL  STRUCTURE  AND  PLATE  TECTONICS  NEAR  THE 
WESTERN  PART  CF  THE  ALASKA  RANGE 

Contribution  by  John  Davies  ar.d  Eduard  Berg 

As  part  of  a  continuing  program  of  crustal  studies  in  Alaska,  5 
temporary,  VHF  telemetered,  seismic  stations  were  operated  at  the  base 
of  the  Alaska  Range  in  the  vicinity  of  Skwentna  (about  120  km  NW  of 
Anchorage),  Alaska,  during  June  and  July  1971.  All  of  the  stations 
were  installed  and  maintained  using  USAF  helicopters  from  Elmendorf, 

Alaska.  The  station  sites  were  chosen  to  increase  the  precision  with 
which  earthquake  hypocenters  could  be  located  in  the  subduction  zone  NW 
of  Anchorage  along  the  Alaska  Range.  The  stations  augmented  those  operated 
routinely  by  the  Geophysical  Institute,  University  of  Alaska,  and  the 

Palmer  Seismological  Observatory,  NOAA. 

Data  from  these  three  sources  is  now  being  analyzed  and  preliminary 

results  promise  that  this  project  will  succeed  in  refining  previous  work 
on  dipping  crust-mantle  interface  and  tectonics  of  the  Cook  Inlet-Susitna 
River  region,  and  be  the  basis  for  further  seismological  studies  in  the 
area.  It  is  important  to  understand  the  seismicity  of  this  region  since 
half  the  population  of  the  State  of  Alaska  lives  in  it. 

This  project  was  significant  because  it  marked  the  first  extensive 
field  use  of  VHF  telemetry  in  seismic  recording  by  our  group.  The  ex¬ 
perience  gained  is  now  being  used  in  other  areas  of  the  State  as  we 
are  operating  about  10  VHF  telemetered  seismic  stations.  The  advantages 
of  this  kind  of  system  are  many  in  a  region  where  stations  are  needed  in 
remote  locations  without  roads,  telephone  channels  or  power. 
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CRUSTAL  TILT  FIELDS  AND  PROPAGATION  VELOCITIES 
ASSOCIATED  WITH  EARTHQUAKES 

A 

Contribution  by  Edvard  Berg  and  William  Lutschak 

ABSTRACT 

Tilt  steps  similar  to  strain  steps  have  been  observed  for  earth¬ 
quakes  with  magnitudes  ranging  from  2  to  8  and  distances  ranging  any¬ 
where  from  10  km  to  12,000  kilometers.  Tilt-step  propagation  velocities 
from  the  hypocenter  to  the  station  have  been  observed  from  1.2  to  3.8 
km/sec,  and  those  corresponding  to  the  arrival  time  of  teleseismic  S 
waves.  Tilt  directions,  amplitudes  and  velocities  observed  at  several 
stations  simultaneously  for  the  same  earthquake  are  internally  consistent 
and  are  likely  to  depend  on  the  tectonic  environment  of  the  station  in 
addition  to  the  focal  mechanism.  The  low  velocity  of  1.2  km/sec  is  con¬ 
sistent  with  a  plastic  wave  propagation.  Other  velocities  are  close  to 
the  Rayleigh  wave  short-period  group  velocity  over  continental  paths. 


* 

Now:  Hawaiian  Institute  of  Geophysics,  University  of  Hawaii 
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INTRODUCTION 


Tilt  and  strain  steps  associated  with  earthquakes  have  been  observed 
widely.  Their  amplitudes  often  exceed  those  theoretically  expected  from 
elastic  fields  associated  with  a  dislocation  at  the  source.  Observations 
on  propagation  velocities  are  sparse  and  mostly  confined  to  a  single  path 
(Wideman  and  Major,  1967;  Berg  and  Pulpan,  1971).  Tilt  (and  strain)  step 
amplitudes  and  velocities  are  presented  for  teleseisms  traveling  through 
different  regions  at  near  S-  and  Rayleigh-wave  velocities,  and  for  local 
earthquakes  recorded  at  several  stations  simultaneously.  Observations  of 
some  local  earthquakes  cover  a  full  quadrant  and  should  provide  some  under¬ 
standing  of  near-source  dynamics  in  relation  to  source  and  tectonic  setting. 

INSTRUMENTATION  AND  METHOD 

The  University  of  Alaska's  Tripartite  Tilt  network  at  MCK,  PAX  and 
GLM  has  been  presented  by  Berg  and  Pulpan  (1971).  The  network  is  com¬ 
prised  of  borehole  installations  of  horizontal  Lamout-Lunar  type  15-second 
pendulums  equipped  with  capacitor  transducers.  The  HIG  (21°18'05"N, 
157°49'20"W)  installation  consists  of  NS  and  EW  Sprengnether  LP  pendulums 
(at  Tq  =  20  seconds)  mounted  on  a  concrete  pier  in  the  basement  of  the 
institute  building  (Hawaii  Institute  of  Geophysics).  Output  at  the  time 
of  the  earthquakes  reported  here  was  recorded  to  DC  from  a  Sprengnether 
capacitor  transducer  on  Helicorder  paper.  The  ALPA  installation  north 
of  Fairbanks  consists  of  triaxial  long-period  (Tq  =  20  seconds)  borehole 
seismometers  with  velocity  transducers  and  position  monitor.  Station 
locations  are  shown  in  figure  15.  Only  a  limited  amount  of  data  has  become 
available  to  far. 

To  determine  amplitudes  and  velocities  for  teleseisms  and  regional 


events,  either  the  digital  records  (1  data  point/second)  or  the  tidal 
outputs  (both  to  DC)  were  used  together  with  the  epicenter  and  origin 
time  of  the  U.S.C.G.S.  preliminary  or  monthly  determinations.  Great 
circle  path  distances  were  determined  through  a  geodetic  computer  pro¬ 
gram  by  Gary  Schnelzer  (Hawaii  Institute  of  Geophysics).  For  the  last 
two  events  in  Table  1,  arrival  times  were  read  from  the  records  presented 
in  the  literature  by  Stacey  and  Rynn  (1970)  and  Tocher  and  Brown  (1971), 
and  their  tilt  and  strain  steps  and  distances  were  used.  (Those  of  Wide- 
man  and  Major,  1967,  were  not.) 

For  local  Fairbanks  earthquakes,  Berg  and  Pulpan  (1971)  presented 
a  number  of  offset  tidal  records  and  the  corresponding  pulses  on  the  LP 
outputs  (including  the  10  October  1970  earthquake  for  GLM) .  GLM  is  lo¬ 
cated  inside  ALPA.  An  analog  output  for  a  step  input  (Cal  pulse)  of 
one  of  the  ALPA  seismometers  was  obtained  (lower  left  corner  of  Figure 
15).  If  the  records  of  local  earthquakes  from  GLM  showed  tilts,  and 
the  pulses  propagating  through  ALPA  could  be  superimposed  on  the  cali¬ 
bration  pulse  form,  the  corresponding  output  was  considered  to  result 
from  a  tilt. 

During  the  9th  International  Symposium  on  Geophysical  Theory  on 
Computers  the  question  was  raised  whether  or  not  such  a  pulse  output 
could  result  from  the  saturation  of  the  ALPA  seismic  amplifiers  by  large 
amplitude  short-period  waves  of  the  local  earthquakes.  However,  this 
seems  unlikely.  In  several  cases  the  tilt  only  occurred  in  one  direc¬ 
tion,  whereas  the  other  perpendicular  output  shows  only  noise.  As  an 
example,  consider  the  NS  component  of  station  34  (Figure  14c)  that  only 
showed  noise  (trace  343).  The  output  level  corresponds  to  3.73  x  10  = 
37.3  mp  at  the  maximum  system  response  (at  25  sec  period,  Table  2, 
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Seismogram  114,  Channel  343).  Since  this  output  is  vectorially  composed 
of  the  three  outputs  of  the  triaxlal  system  it  should  show  saturation  at 
a  much  higher  level  (comparable  to  the  pulse  amplitude  of  the  east-west 
output  (Figure  14b)  trace  342)  if  saturation  of  one  or  more  amplifiers 
had  taken  place.  Since  it  does  not  show  such  a  saturation  it  is  con¬ 
cluded  that  the  pulse  shown  on  the  east-west  output  is  real  (Figure  14b, 
trace  342)  corresponding  to  an  output  level  of  2.2  y  at  the  maximum  system 
response  at  25  sec  period  (Table  2,  seismogram  113,  Channel  342).  Simi¬ 
larly  it  is  felt  that  the  other  pulses  that  are  like  a  step  input  to  the 
calibration  coil  (shown  in  the  lower  left  of  Figure  15)  are  real.  In 
other  cases  the  pulse  is  clearly  visible  but  of  smaller  amplitude  than 
the  maximum  noise  level  recorded  during  the  time  interval  for  which  the 
record  has  been  obtained  (as  station  36  for  12  November  1970,  Figure  19). 

For  some  ALFA  outputs  with  responses  deviating  somewhat  from  the  ideal 
step  response,  overshoot  amplitudes  were  used  to  determine  tilt  ampli¬ 
tudes.  Absolute  tilt  amplitudes  were  obtained  by  converting  the  given 
output  amplitudes  via  calibration  amplitudes  (or  voltages)  and  using 
mgsina  %  mga  =  iG,  where  m  =  mass  of  seismometer,  g  earth  gravity,  a  = 
deviation  angle  of  seismometer  boom  for  the  equivalent  pendulum  length 
of  a  20  sec  pendulum  when  the  calibration  current  i  is  applied  and  G  is 
the  calibration  coil  constant.  Time  resolution  on  the  computer-plots 
was  1  inch/100  seconds  for  the  10  October  1970  and  12  November  1970 
earthquakes  and  2  inches/100  seconds  for  the  others.  Original  computer 
outputs  are  presented  in  Figures  14  and  15  near  corresponding  station 
locations  and  are  so  oriented  that  ground  motion  N  and  E  corresponds 
to  record  trace  N  and  E.  The  corresponding  tilt  amplitude  for  the  trace 
(normalized  to  +  1  inch)  is  also  shown.  The  resulting  tilt  direction 
and  amplitude  are  presented  in  the  figures.  Results  for  COL  have  been 


similarly  obtained. 

Only  those  records  that  matched  the  calibration  pulse  have  been  used 
to  determine  the  arrival  times  shown  in  Figure  10  for  local  quakes  (ex¬ 
ceptions  or  poor  results  for  regional  and  distant  earthquakes  are  given 
in  brackets  in  Table  1) . 

Focal  point  and  origin  time  have  been  determined  from  a  local  short- 
period  telemeter  network  using  the  S-P  t.  mes  at  LV  and  P  arrivals  at  all 
other  short-period  stations  (see  Figures  15  to  19  for  locations) .  S-P 
times  at  LV  for  the  four  earthquakes  used  in  Figure  5  ranged  from  1.8  to 
3.7  seconds. 


RESULTS 

A)  Distant  and  Regional  Earthquakes 

Tilt  (and  strain)  amplitudes  and  propagation  velocities  for  a  distance 
range  of  70  km  to  12,600  km  and  a  magnitude  range  of  3.6  to  8.0  are  pre¬ 
sented  in  Table  1.  The  velocities  determined  from  arrivals  on  digital  or 
Helicorder  records  should  be  considered  very  reliable,  whereas  those  ob¬ 
tained  from  tidal  records  (25  mm/hour)  are  presented  only  as  being  indica¬ 
tive.  Three  rather  distinct  velocity  groupings  are  apparent:  that  of  2.6 
to  2.8  km/sec  traveling  through  parts  of  the  Aleutian  arc  (6  November  1971), 
South  Central  Alaska  (26  March,  1971),  and  in  an  area  north  of  Fairbanks 
(11  May  1971);  that  near  3.1  km/sec  traveling  at  least  partially  through 
the  continental  area;  and  that  of  3.5  to  3.8  km/sec  traveling  through  the 
Pacific.  All  three  types  of  step  arrivals  seem  to  be  associated  with  and 
somewhat  follow,  in  time,  the  maximum  amplitudes  of  the  Rayleigh  waves. 

A  fourth  step  arrival  follows  closely  the  S-wave  arrival:  it  is  visible 
on  the  Stacey  and  Rynn  (1970,  their  Figure  1)  illustration  (17  January 
1968)  and  on  our  digital  and  Helicorder  records  (10  January  1970,  29 
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February  1972).  A  tilt  associated  with  the  P-wave  arrival  might  also  be 

present  in  the  Stacey  and  Rynn  record. 

It  should  be  noted  that  tilt  steps  at  PAX  and  MCK.  for  the  New  Guinea 
and  the  Solomon  Islands  earthquakes  <10  January  1971  and  14  July  1971, 
respectively)  are  nearly  identical  in  direction  and  amplitude  (Figures  5 
and  6).  The  stations  are  nearly  200  km  apart,  along  the  Denali  Fault  in 
the  Alaska  Range,  whereas  at  GLM,  on  the  northern  edge  of  the  Tanana  Basin 
250  km  away,  the  tilt  steps  show  much  smaller  amplitudes,  for  roughly  simi¬ 
lar  epicentral  distance  and  azimuth.  Both  earthquakes  also  had  comparable 

magnitudes. 

The  rise  time  of  the  tilt  steps  seems  to  be  small  when  compared  to 
seismometer  and  circuit  time  constants.  Most  of  the  tilt  steps  do  not  seem 
to  recover  with  time.  Figures  1  through  9  show  typical  tidal  and  digital 
records . 

B)  Local  Earthquakes 

Results  of  tilt-step  arrival  times  versus  focal  distance  are  presented 

in  Figure  10.  The  earthquake  magnitudes  ranged  from  3.0  upward  and  depths 

were  18  to  20  km  and  11  km  (1  May  1971).  Other  observations  not  presented 

here  include  an  earthquake  of  magnitude  2.4  and  depth  of  14  km,  for  which 

* 

tilt  steps  were  recorded  at  GLM  and  only  at  nearby  ALPA  stations.  Since 
the  stations  are  very  near  the  epicenters,  focal  distance-travel  time 
served  to  calculate  the  step  velocities, 

*A  magnitude  2  earthquake  was  the  smallest,  for  which  a  tilt  step  (  of  2 
msec  arc)  has  been  observed  at  GLM. 
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Three  groups  of  velocities  are  indicated:  one  of  aVout  4  km/sec  which 
includes  the  3.6  km/sec  (10  October  1970,  because  of  relatively  poor  time 
resolution  on  the  computer  plots);  a  second  of  about  2.3  to  2.5  km/sec;  and 
the  third  of  1.2  km/sec  (13  May  1971).  It  should  be  noted  that  for  the 
earthquake  of  1  May  1971,  stations  33,  A2  and  32  are  all  aligned  with  the 
epicenter  and  GLM  is  only  slightly  off.  If  the  2.40  km/sec  line  through 
the  arrivals  at  those  stations  is  adjusted  i'o  as  to  pass  0  focal  distance 
at  a  time  different  from  the  origin  time,  a  better  fit  is  obtained,  resulting 
in  a  somewhat  lower  velocity  near  2.2  km/sec. 

Rise  times  on  the  ALPA  output  correspond  (from  visual  matching  by  super¬ 
position)  to  that  of  a  calibration  step  input,  and  at  GLM  are  small  compared 
to  instrument  and  filter-time  constants.  For  GLM  ;his  has  also  been  verified 
for  unfiltered  digital  outputs  from  local  earthquakes  (Figure  1). 

C)  Amplitude-Distance 

The  combined  vectorial  amplitude  of  the  tilt  steps  for  the  ALPA,  GLM 
and  COL  stations  are  presented  for  7  earthquakes  that  originate  near  these 
stations.  Figures  11a,  b  show  tilt  amplitudes  versus  epicentral  distance; 
all  data  are  summarized  in  Figure  12.  Tilt  amplitude  and  direction  of  5  of 
the  local  quakes  are  presented  in  Figures  15  through  19.  The  numbers  at 
the  tip  of  the  tilt  arrow  indicate  the  factor  by  which  the  length  of  the 
arrow  has  to  be  multiplied  to  obtain  tilt  amplitude  corresponding  to  the 
scale  in  the  right  hand  lower  corner  of  the  figures.  Local  magnitudes 
range  from  2.8  to  somewhat  over  four  and  have  been  determined  from  short- 
period  vertical  station  LV,  BLR  and  SCM.  LV  is  located  in  the  epicentral 
area,  BLR  is  some  150  to  180  km  and  SCM  over  300  km  from  the  epicenters. 
Epicentral  distance  for  the  tilt  steps  recorded  range  from  8  to  80  km. 

Data  points  (tilt  amplitude  versus  epicenter  distance)  in  brackets  in¬ 
dicate  either  possible  maximum  values  (especially  below  10  msec  arc)  or 
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are  doubtful.  Stations  indicated  by  the  same  subsymbol  (vertical,  hori¬ 
zontal  or  oblique  line)  are  located  on  or  nearly  on  a  line  passing  through 
the  epicenter  indicated  by  a  circle  and  cross  in  Figures  15  through  19. 
Amplitude  die-off  with  distance  R  as  R-X  is  as  indicated  next  to  the  lines 
joining  corresponding  in  line  stations.  The  value  for  x  varies  from  1.6 
to  over  7.  For  none  of  the  7  earthquakes  all  ALFA  data  have  been  obtained. 
However,  in  many  instances  stations  did  not  show  any  pulse  or  other  motion 
related  to  the  earthquakes,  out  of  the  noise  level.  The  noise  level  indi¬ 
cated  by  the  output  was  lowest  for  the  three  northern  most  stations.  Those 
stations  that  showed  only  noise  are  indicated  in  Figure  19  by  an  additional 
circle  around  the  station  location.  There  is  doubt  concerning  the  output 
of  station  11  for  the  May  1971  quakes,  showing  exceedingly  high  amplitude, 
compared  to  surrounding  stations  or  earlier  earthquakes.  The  high  amplitude 
is  due  to  a  single  component  of  the  triaxial  system  (aligned  with  the  struc¬ 
tural  trends  in  the  area) ,  whereas  the  two  other  components  show  very  small 
amplitudes  on  the  raw  data.  A  similar  sftuation  occurs,  however,  for  other 
ALPA  stations  showing  "reasonable"  outputs,  a  fact  that  would  argue  for 
retention  of  data  for  the  May  1971  quakes  at  station  11. 


DISCUSSION 


A)  Distant  and  Regional  Earthquakes 

I 

Maue  (1954)  has  investigated  (two  dimensionally)  the  motion  of  re- 

I 

laxauion  that  occurs  when  an  elastic  medium  under  tension  (in  y  direction) 
is  suddenly  cut  open  from  one  side,  perpendicular  to  the  tension  axis 
(along  x  <  0,  y  =  0) .  The  plane  compressional  wave  corresponding  to  a 


step  is  the  strain  at  the  shockfront,  and  the  amplitude  is  such  that 
the  initial  tension  is  compensated.  On  the  other  hand,  if  a  solution 
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corresponding  to  a  sinusoidal  point  force  in  a  semi-infinite  medium  may 
be  expressed  as  g  (m,  coordinates,  elastic  and  source  parameters)  e*0^  , 
where  t*  is  the  time  corresponding  to  a  particular  propagation  mode,  and 
a  step  in  force  is  applied  at  the  source  at  t  =  0,  the  resulting  solution 
is  given  by  the  fourier  integral  over  the  fourier  transform  of  the  source 
time  function  times  the  solution  function  (Ewing  et  al.,  1957,  page  61) 
integrated  over  w.  Therefore  a  step  in  stress  (and  strain  and  tilts)  is 
expected  at  the  recording  site  at  times  corresponding  to  the  particular 
propagation  modes.  The  observations  of  steps  traveling  at  velocities  near 
S-  and  Rayeligh-wave  velocities  seem  to  confirm  this  contention. 

Moreover,  the  velocities  of  steps  occurring  during  the  Rayleigh-wave 
train  seem  to  correspond  closely  to  the  high-frequency  portion  of  the  con¬ 
tinental  Rayleigh  waves  group  velocity  (i.e.,  3.1  km/sec),  and  for  the 
local  earthquakes,  also  close  to  the  velocities  of  sedimentary  Love  and 
Rayleigh  waves  (2.0  -  2.5  km/sec)  (Oliver,  1962).  The  velocities  (Table 
1)  across  parts  of  the  Pacific  closely  correspond  to  those  observed  by 
Wideman  and  Major  (19u7)  for  oceanic  paths. 

The  velocities  of  tilt  (  and  strain)  steps  associated  with  the  S  waves 
correspond  closely  to  the  S-wave  arrival  time.  Particularly  noteworthy 
is  the  observation  of  a  compressional  strain  in  the  NW-SE  direction  asso¬ 
ciated  with  such  an  S-wave  arrival  time  from  the  Amchitka  explosion  at  Cape 
Sarichef  (record  in  Tocher  and  Brown,  1972).  Calculations  of  the  strain 
and  tilt  field  resulting  from  a  purely  elastic  dislocation  model  at  the 
source  have  been  carried  out  by  different  authors.  However,  the  tilt  am¬ 
plitudes  for  the  teleseismic  events  often  exceed  the  theoretically  expected 
ones  by  many  times  (as  in  Press,  1965)  .  Moreover,  there  also  are  large 
amplitude  differences  among  PAX  and  MCK  on  the  one  hand  and  GIM  on  the 
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other  for  the  New  Guinea  and  the  Solomon  Island  earthquakes  for  similar 
epicentral  distance  and  azimuth.  For  a  different  azimuth  (the  Chile  earth¬ 
quake),  however,  amplitude  differences  (for  comparable  magnitudes)  are  much 
smaller  at  PAX,  MCK  and  their  relation  to  those  at  GLM  is  different.  This 
supports  Stacey  and  Rynn's  idea  that  the  tilt  steps  are  related  to  local 
tectonics  -  in  the  Alaska  case  on  a  large  scale,  since  the  records  for  MCK 
and  PAX  near  the  Alaska  Range  are  remarkably  similar,  whereas  those  at  GLM 
north  of  the  Tanana  Basin  are  different  (compare  Figures  5  through  9).  It 
is  too  early  to  compare  the  tilt-step  direction  to  a  local  secular  trend, 
because  recording  times  are  not  long  enough.  A  more  detailed  discussion 
of  the  amplitude  fall-off  with  distance  has  been  given  in  an  earlier  paper 
(Berg  and  Pulpan,  1971) . 

B)  Local  Earthquakes 

Any  discussion  of  the  results  of  local  earthquakes  is  a  most  ambiguous 
one  since  only  one  full  quadrant  of  data  for  the  deformation  of  the  crustal 
material  in  the  form  of  tilt  steps  and  their  arrival  times  is  available,  and 
theoretical  calculations  for  surface  deformation  resulting  from  strike  slip 
faulting  at  distances,  comparable  to  source  depth,  and  realistic  velocity- 
depth  functions  are  almost  nonexistent. 

(a)  Tilt  Direction 

The  assumption  of  almost  vertical  faults  with  strike  slip  motion  re¬ 
sulting  from  a  N-  to  NW-oriented,  compressional  axis  in  the  area  is  sup¬ 
ported  by  short-period  and  tilt  investigations  (Berg  and  Pulpan,  1971; 

Gedney  and  Berg,  1969a  and  1969b).  The  tilt  field  directions  of  the  10 
October  1970  earthquake  (Figure  18)  are  consistent  with  theoretical  cal¬ 
culations  (Chinnery,  1961;  Press,  1965)  for  a  serai-infinite  medium  if  the 
line  for  zero  vertical  uplift  (that  is  roughly  at  a  distance  from  the 
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epicenter  which  is  cluse  to  the  focal  depth)  is  placed  beyond  stations 
BRH  and  33.  In  this  case,  all  short-period  first-motion  observations 
(Berg  and  Pulpan,  1971)  at  MCK,  HPP,  LV,  MIN,  BRH,  GLM  and  MCB  are  con¬ 
sistent,  PJD  actually  showing  only  a  small  up-motion  followed  by  a  large 
down-motion,  indicative  of  a  location  near  the  nodal  plane.  Of  the  ALPA 
stations,  only  station  3A  shows  a  downward  motion  inconsistent  with  the 
interpretation,  whereas  the  down-motion  at  station  23  and  the  up-motion 
at  station  33  might  result  from  nearness  to  the  nodal  line  and/or  a  small 
dip  in  the  nodal  plane.  The  nodal  planes  required  to  fit  this  solution 
are  nearly  NS  and  EW,  the  NS  plane  passing  slightly  to  the  east  of  GLM, 
and  the  motion  is  either  no,;th  on  the  eastern  side  of  the  NS  plane  or  east 
on  the  northern  side  of  the  EW  plane,  consistent  with  a  northeasterly  com- 
pressional  axis.  It  is  believed  that  tilt-field  directions  would  not  be 
altered  drastically  if  a  realistic  velocity-depth  function  were  to  be  used 
in  the  theoretical  calculations,  and  the  amplitudes  might  be  closer  to 
observations.  R.  Sato  (personal  communication,  April,  1972)  has  shown 
theoretically  that  - ae  Inclusion  of  a  1.7  km  thick  surface  layer  (V  ■  3.5ft, 

Vs  "  2,37  km/ sec)  followed  by  a  2.8  km  thick  layer  (V  *  5.5,  V  »  3.10 

P  s 

km/sec)  over  higher  velocity  layers  in  the  crust  drastically  alters  the 
deformation  surface  field  fot  dip-slip  and  inclined  strike-slip  dislo¬ 
cations,  whereas  numerical  calculation  for  a  model  with  a  first  layer  of 
33  km  thickness  (Vp  «*  6.3,  Vg  ®  3.55  km/sec)  results  in  surface  deformation 
closely  similar  to  that  of  a  semi-infinite  medium  (Sato,  1971). 

(b)  Tilt  Velocities 

The  arrival  times  most  often  observed  for  focal  distances  to  40  km 
correspond  to  a  velocity  of  2.0  to  2.3  km/sec.  These  might  be  interpreted 
as  associated  with  the  Rayleigh  wave  for  a  sedimentary  layer  at  the  edge 
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Of  the  Tanana  Basin.  Hanson  at  al.  (1968)  found  such  a  layer  2.1  to  2.6 
km  thick  (Vp  -  3.67,  Vs  ■  2.31  km/sec)  at  least  to  the  edge  of  the  basin. 
However,  parts  of  the  travel  path  are  through  Birch  Creek  schist  and  the 
Interpretation  is  open  to  question.  At  larger  distances  (beyond  70  km) 
velocities  near  2.8  to  3.1  km/sec  have  been  observed  corresponding  to  those 
for  the  pulses  associated  with  Rayleigh  waves  over  continental  paths  (see 
Figure  10,  13  May  1971,  station  31). 

No  explanation  is  offered  for  the  4  km/sec  branch  which  shows  a  delay 
of  nearly  10  seconds  for  the  10  October  1970  and  16  January  1971  earthquakes 
but  no  delay  for  the  shallow-depth  earthquake  of  1  May  1971.  A  very  low 
propagation  velocity  (1.2  km/sec)  was  observed  between  stations  33  and  A2 
(in  line  with  the  epicenter  of  the  13  Hay  1971  earthquake)  and  occasionally 
for  arrivals  at  a  single  station  (like  that  at  station  24,  16  January  1971). 
The  lowest  velocity  observed  was  0.9  km/sec.  Likewise  Romig  and  others 
(1969,  their  Figure  3.  site  2)  seems  to  Indicate  a  strain  step  arrival  about 
1  minute  to  1  minute  and  20  sec  after  the  origin  time  of  BENHAM  ,t  a  dis¬ 
tance  of  71.5  km.  Indicating  a  similar  low  velocity  of  propagation  near 
1  km/sec.  These  low  velocities  eecm  to  be  associated  with  the  propagation 
of  a  plastic  pulse  l„  a  highly  stressed  area.  By  superposing  the  earth¬ 
quake-generated  stress,  the  nonlinear  portion  of  the  stress-strain  curve 

is  responsible  for  a  plastic  wave  that  will  travel  with  a  velocity  V  as  given 
by  Kolsky  (1953): 

y2  '  <  3f  /  o  ) 

where  e  -  stress,  c  -  strain  and  p  -  density.  In  the  linear  part  of  the 
.trees  strain  curve  d,/dc  -  E  or  Young’s  modulus.  In  the  nonlinear  part, 
for  V  -  1.2  km/sec  and  p  .  2.7  gm/cm3,  a  value  of  3.9  x  1010  dyne/cm2  is 
obtained  for  the  modulus.  This  i,  twenty  times  lower  than  the  corresponding 
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E  modulus  for  a  material  with  Vp  •=  5.5  km/sec.  As  a  consequence,  such  a 
plastic  wave  will  generate  a  considerably  larger  deformation  than  the 
elastic  wave  for  the  same  applied  stress.  A  plastic  deformation  could 
also  be  triggered  by  stresses  from  teleseismic  sources  in  a  prestressed 
area  and  could  perhaps  alter  the  direction  of  the  tilt  (and  strain  axis) 
and  account  for  some  of  the  larger  tilt-amplitudes  observed  in  tectonically 
active  areas,  such  as  the  Alaska  Range. 

(c)  Tilt  Amplitudes 

A  number  of  features  emerge  from  the  data.  Most  of  the  local  earth¬ 
quakes  are  at  a  depth  near  20  km,  so  that  the  theoretically  strong  depen¬ 
dence  of  deformation  amplitudes  on  depth  (like  H_2)  is  not  considered  in 
this  discussion. 

Tilt  amplitudes  for  the  same  distance  and  magnitude  scatter  at  least 
by  an  order  of  magnitude  and  seem  dependent  on  station  location.  As  an 
example,  station  24  shows  about  one  order  less  in  amplitude  (compared  to 
other  ALPA  stations)  for  similar  epicentral  distance  (Figure  12) .  Station 
31  seems  to  indicate  higher  amplitudes  (when  recorded)  that  might  be  struc¬ 
turally  controlled,  indicated  by  the  predominantly  WSW  tilt  direction  for 
the  Fairbanks  quakes  and  almost  opposite  direction  for  the  Min to  quake  at 
the  western  side  of  ALPA  (see  Figures  15  through  19) . 

A  general  increase  in  amplitude  with  magnitude  is  indicated  at  least 
to  a  distance  between  30  to  50  km  from  the  epicenter.  This  is  demonstrated 
by  Figure  13  for  stations  34,  24  and  COL,  where  tilt  amplitude  versus 
magnitude  is  plotted.  These  particular  stations  are  located  in  a  narrow 
epicentral  distance  range  (35.6  to  39  km  for  24,  32.8  to  38.2  km  for  34, 

12  and  15  km  for  COL)  except  for  the  10  Oct  70  mag.  3.2  quake  (where  24 
and  34  are  near  51  km).  The  two  shallower  earthquakes  indicate  somewhat 
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higher  amplitude  for  24  (two  data  points  above  full  line,  Figure  13)  but 
not  for  34,  where  the  amplitude  for  the  mag  2.9  quake  is  very  low.  For 
the  dashed  lines  the  log  of  amplitudes  shows  direct  dependence  on  M, 
whereas  for  the  full  lines  dependence  is  like  M/1.5,  a  result  in  line 
with  earlier  findings  by  Berg  and  Pulpan  (1970) . 

A  cutoff  magnitude  for  recording  tilt  steps  seems  at  least  indicated 
for  34  (at  a  distance  near  35  km) .  Taking  all  data  (of  Figure  12)  together 
it  is  not  clearly  possible  to  ascertain  a  dependence  of  a  cutoff  distance 
on  Magnitude.  Some  cutoff  between  40  and  50  km  is,  however,  indicated  by 
the  reliable  data,  roughly  at  a  distance  twice  the  source  depth. 

CONCLUSION 

Observed  tilt  steps  often  seem  to  be  too  large  to  be  accommodated  by 
existing  theoretical  dislocation  models  for  a  single  layered  homogeneous 
isotropic  and  elastic  crust.  Alternative  solutions  should  be  looked  for 
in  terms  of  more  realistic  crustal  models  and  include  plastic  deformation 
in  tectonically  stressed  areas,  related  to  the  nonlinear  portion  of  the 
stress  strain  curve.  This  is  partially  indicated  by  the  large  tilt  am- 
plitude  at  short  distances  from  the  epicenter  and  by  the  observed  low  pro¬ 
pagation  velocity  of  tilt  steps  for  local  earthquakes  in  the  Fairbanks  area, 
and  partially  by  the  high  tilt  amplitudes  associated  with  large  distant 
earthquakes  in  the  tectonically  active  area  of  Alaska. 
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Figure  1. 
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GLM  LP  digital  records  from  Min to  earthquake  unfiltered; 
calibration  pulse  box  car  filtered  at  5  seconds.  Up  on 
records  is:  Groundmotion  W  and  N;  Tilt  E  and  S.  Original 
scale  was  3/4  inch  per  minute. 

PAX  digital  records  from  an  earthquake  in  the  Gulf  of  Alaska. 
Left  side  unfiltered  data,  right  side  sensitivity  increased 
five  times  and  data  are  box  car  filtered  at  20  sec.  Notice 
tilt  offset  towards  the  southeast.  Time  scale  on  original 
records  is  3/4  inch  per  minute.  Indicated  time  is  not  correct. 

Unfiltered  digital  data,  PAX,  for  the  New  Guinea  quake  of 
10  January,  .1971,  distance  over  9000  km,  original  time  scale 
1/4  inch  per  minute.  Compare,  to  next  figure. 

,  PAX  record  of  New  Guinea  earthquake  10  January  1971  box  car 
filtered  at  100  sec  and  5  times  increased  sensitivity  grom 
previous  figure.  Notice  the  firs£  tilt  offset  near  07  40 
toward  N-E  and  the  second  near  08  08.  Propagation  velocity 
the  second  tilt  is  3.15  +  0.1  km/sec.  Compare  also  to  next 
figure . 

.  E-W  tilt  record  of  the  January  10,  1971  Guinea  earthquake, 
recorded  at  MCK,  GLM  and  PAX.  The  difference  in  amplitudes 
results  from  the  difference  in  operation:  MCK  was  directly 
recorded  at  the  station,  GLM  was  telemetered  and  (R-C)  fil¬ 
tered  at  a  30  sec  corner  period.  Both  traces  are  feedback 
signals.  PAX  is  the  Transducer  output  (no  feedback  applied), 
telemetered  and  recorded  after  a  (R-C)  30  sec  corner  period 
filter.  "0"  reference  marks  are  spaced  1  hour  apart.  Up  on 
record  is  tilt  toward  the  indicated  direction,  and  the  bars 
correspond  to  50  msec  arc. 

.  14  July  1971  Solomon  Island  Earthquake.  E-W  tilt  records. 

Compare  the  very  similar  tilts  during  the  New  Guinea  earth¬ 
quake  in  the  preceeding  figure.  Note  however  that  recording 
polarity  for  MCK  is  reversed,  but  tilt  direction  is  not.  Other 
details  as  in  Figure  5. 

.  9  July  1971  Chile  Earthquake,  E-W  tilt  records.  For  other 

details  refer  to  Figure  5. 

.  14  July  1971  Solomon  Island  N-S  tilt  records.  For  other  de¬ 

tails  refer  to  Figure  5. 

.  9  July  1971  Chile  Earthquake  N-S  tilt  records.  For  other 

details  refer  to  Figure  5. 


Figure  10. 

Figure  11. 

Figure  12. 
Figure  13. 

Figure  14. 

Figure  14. 

Figure  14. 
Figure  15. 

Figure  16. 

Figure  17. 

Figure  18. 

Figure  19. 


Tilt  arrival  times  after  origin  time  versus  focal  distance 
for  four  Fairbanks  earthquakes.  Station  identification  g  ven 

by  numbers  and/or  letters. 

a.)  b.)  Total  tilt  amplitude  versus  epicentral  distance  values 
in  (  )  are  uncertain  or  possible  maximum.  Stations  with  t  e 

same  subsymbol  (horizontal  or  vertical  lines)  are  located  in 
or  nearly  in  line  with  the  epicenter.  Plots  are  ordered  by 
decreasing  magnitude. 


Summary  plot  of  previous  Figure  11a  and  b. 

Tilt  amplitude  versus  magnitude  for  stations  34,  24  and  COL; 
each  station  is  at  a  comparable  epicenter  distance  for  all 
quakes,  and  focal  distance  was  also  nearly  equal  except  for 
2  quakes  (see  text) . 

a. )  Vertical  ALPA  outputs  for  10  October  1970  quake.  The  first 
two  letters  and/or  numbers  are  station  identification. 

b. )  E-W  ALPA  outputs  for  10  October  1970  quake  leading  to  tilts 
indicated  in  Figure  19.  For  amplitudes  refer  to  Table  2. 


c.)  N-S  ALPA  outputs  for  10  October  1970  quake. 

Normalized  ALPA  station  outputs,  tilt  directions  and  tilt 
amplitudes  (including  GLM  and  COL)  for  the  1  May  1971  quake. 

For  comparison  with  calibration  pulse  refer  to  text.  The 
length  of  tilt  arrow  has  to  be  multiplied  by  the  factor  in 
dicated  to  correspond  to  scale. 

Tilt  direction  and  amplitude  for  the  11  May  1971  quake.  The 
length  of  the  arrow  has  to  be  multiplied  by  the  factor  indicated 

to  correspond  to  scale. 

Tilt  direction  and  amplitudes  for  the  13  May  1971  quake.  The 
length  of  the  arrow  has  to  ba  multiplied  by  the  factor  indicated 

to  correspond  to  scale. 

Tilt  direction  and  amplitudes  for  the  10  October  1970  quake. 

The  length  of  the  arrow  has  to  be  multiplied  by  the  factor 
indicated  to  correspond  to  scale.  The  original  data  are  pre 
sented  in  Figure  14  and  amplitudes  are  given  in  Table  2. 


Tilt  direction  and  amplitude  for  the  12  November  1970  quake. 
The  length  of  the  arrow  has  to  be  multiplied  by  the  factor 
indicated  to  correspond  to  scale. 
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CHANNEL  10  MAT  11  1971  NO  FILTEfl  £00  .4V/INCH  M.INTO  68  KM  MAG  36 


Glm  LP  digital  records  from  Ilinto  earthquake  unfiltered; 
Calibration  pulse  box  car  filtered  at  5  seconds.  Up  on 
records  is:  Groundnotion  W  and  11;  Tilt  E  and  S.  Original 


CHANNEL  16  MATCH  zl  1971  SDO  MV/ INCH 


NO  FRTEM 


FiR.  2  PAX  digital  records  from  an  earthquake  in  the  Gulf  of 

Alaska,  Left  side  unfiltered  data,  rip,ht  side,  sensitivity 
is  increased  five  times  and  data  are  box  car  filtered  at 
.  20  sec.  Wotice  tilt  offset  towards  the  southeast.  Time 

scale  on  original  records  is  3/4  inch  per  minute.  Indicated 
time  is  not  correct. 
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E-W  tilt  record  of  the  January  10,  1971  Guinea  earthquake, 
recorded  at  MCK,  GUI,  and  PAX.  The  difference  in  amplitude 
results  from  the  difference  in  operation:  MCK  was  direct  y 
recorded  at  the  station,  GUI  was  telemetered  and 
filtered  at  a  30  sec  corner  period.  Both  traces  are  feedback 
sirnals.  PAX  is  the  transducer  output  (no  feedback  applied), 
telemetered  and  recorded  after  a  (R-C)  30  sec  corner  period 
filter.  "0"  reference  marks  are  spaced  I  .war  apart.  Up  on 
record  is  tilt  toward  the  indicated  direction,  and  the  bares 

correspond  to  50  msec  arc. 
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Fig.  6  July,  1971  Solomon  Island  earthquake.  E-W  tilt  records. 

Compare  the  very  similar  tilts  during  the  New  Guinea  earth" 
quake  in  the  proceeding  figure.  Note  however  that  recording 
.  polarity  for  IICK  is  reversed,  but  tilt  direction  is  not. 

Other  details  as  in  figure  3. 
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9  July,  1971  Chile  earthquake  N-S  tilt  records.  For  other 
details  refer  to  figure  5. 
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Fig.  11  a.)  Total  tilt  amplitude  versus  epicentral  distance  values 

in  ()  are  uncertain  or  possible  maximum.  Stations  with  the 
same  subsymbol  (horizontal  or  vertical  lines)  are  located 
in  or  nearly  in  line  with  the  epicenter.  Plots  are  ordered 
by  decreasing  magnitude. 


b.)  Total  tilt  amplitude  versus  epicentral  distance  value: 
in  ()  are  uncertain  or  possible  maximum.  Stations  with  thi 
same  subsymbol  (horizontal  or  vertical  lines)  are  located 


m  SEC  ARC 
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Fig.  12  Summary  plot  of  previous  Figure  II  a  and  b. 


I  I  I _ I _ I _ I 

2.5  3.0  3.5  4.0  4.5  MAG 

Tilt  amplitude  versus  magnitude  for  stations  34,  24  and  COL. 
Each  station  is  rt  a  comparable  epicenter  distance  for  all 
quakes,  and  focal  distance  was  also  nearly  equal  except  for 
2  quakes  (see  text). 
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Fig.  14  b.J  E-U  ALPA  outputa  for  10  Octobor,  1970  quaJ-.c  leading  to 
tUl*  ItlrJlcatad  In  Fig,  19.  For  aipiUCtlilQl  nfar  to 
labia  2. 
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Fig.  15  Normalized  ALPA  station  outputs,  tilt  directions  and  tilt 

amplitudes  (including  GUI  and  COL)  for  the  I  May,  1971  quake 
For  comparison  with  calibration  pulse  refer  to  text.  The 
length  of  tilt  arrow  has  to  be  multiplied  by  the  factor 
indicated  to  correspond  to  scale. 


Fig.  16  Tilt  direction  and  amplitude  for  the  II  May,  1971  quake.  The 
length  of  the  arrow  has  to  be  multiplied  by  the  factor 
indicated  to  correspond  to  scale. 
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FIr.  18  Tilt  direction  and  amplitudes  for  the  10  October,  1970  quake. 
The  length  of  the  arrow  has  to  be  multiplied  by  the  factor 
indicated  to  correspond  to  scale.  The  original  data  are 
presented  in  Fip.14  and  amplitudes  arc  Riven  in  Table  2. 
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Fig,  19  Tilt  direction  and  amplitude  for  the  12  November,  1970  quake. 

The  length  of  the  arrow  has  to  be  multiplied  by  the  factor 
indicated  to  correspond  to  scale. 
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